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Abstract 

We consider the anisotropic evolution of spatial dimensions and the stabilization of internal 
dimensions in the framework of brane gas cosmology. We observe that the bulk RR field can give 
an effective potential which prevents the internal subvolume from collapsing. For a combination of 
(D — 3)-brane gas wrapping the extra dimensions and 4-form RR flux in the unwrapped dimensions, 
it is possible that the wrapped subvolume has an oscillating solution around the minimum of 
the effective potential while the unwrapped subvolume expands monotonically. The flux gives a 
logarithmic bounce to the effective potential of the internal dimensions. 
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I. INTRODUCTION 



There has been considerable effort during the past decade in the area of higher dimen- 
sional gravity where the four- dimensional gravity can be recovered as an effective theory. In 
the theory of general relativity, the dimensionality of the universe is an assumption. It is 
not derived dynamically from a fundamental theory. However, unifying gravity with other 
forces of nature strongly suggests that there may be more than three spatial dimensions. 
One plausible explanation can be given by string theory, which is higher dimensional. A 
mechanism to generate three spatial dimensions was proposed in Refs. fl and Q. This 
is a mechanism to generate dynamically the spatial dimensionality of spacetime and to ex- 
plain the problem of initial singularity. The key ingredient of this model is based on the 
symmetry of string theory, called T-duality. With this symmetry, spacetime has a topology 
of a nine- dimensional torus, and its dynamics is driven by a gas of fundamental strings. 
Since strings have (l+l)-dimensional world volumes, they can intersect in 2(1+1) spacetime 
dimensions or less from the topological point of view. The winding modes can annihilate 
the anti- winding modes if they meet. Once the windings are removed from the dimensions, 
these dimensions can expand into the large directions that we observe today. Thus, with a 
gas of strings, three spatial dimensions can become large. 

With the realization of D-branes in string theory [3|, there had been attempts to study 
string cosmology with a gas of D-branes including strings . The mechanism of Branden- 
berger and Vafa (BV) in Refs. l|, Q] was reconsidered in this new framework of string theory. 



Authors of Ref. [6| argued that when branes are included, strings still dominate the evolu- 



;ion of the universe at 



ate times. This model of brane gas cosmology was studied extensively 
1^ | (see also 15] for comprehensive reviews). In the picture of brane 



gas cosmology, the universe starts from a hot, dense gas of D-branes in thermal equilibrium. 
The winding modes of branes prevent the spatial dimensions from growing. Branes with 
opposite winding numbers can annihilate if their world volumes intersect. Thus, a hierarchy 
of scales can be achieved between wrapped and unwrapped dimensions. Though the BV 
mechanism can trigger the onset of cosmological evolution based on string theory, this does 
not guarantee the answer to the question why there are three large spatial dimensions while 
others remain small. There have been attempts to yield the anisotropic expansion based on 
certain wrappings of the brane gases 10, 1 1 it ] - 
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It is well known that the stability of internal dimensions cannot be achieved with string 
or brane gases of purely winding modes. Stabilizing the extra dimensions is an important 
issue in string or brane gas cosmology [l2, 14. Im \v\ . For a particular configuration 
of brane gas in eleven dimensional supergravity, Campos [13] has shown the importance of 
fluxes in stabilizing the internal dimensions through numerical computations of the evolution 
equations for wrapped and unwrapped dimensions. It is shown that all internal dimensions, 
except the dilaton, are dynamically stabilized in the presence of two-form flux with a gas of 



long strings 14]. In certain string theories 171 ]. special string states can stabilize the extra 



dimensions for fixed dilaton. Massless string modes at self-dual radius play a crucial role for 
the radion stabilization. 

In this work we will focus on the anisotropic evolution of spatial dimensions and the 
stabilization of internal dimensions based on a certain configuration of branes and bulk 
background field. We consider the possible stabilization of internal dimensions with bulk 
Ramond-Ramond (RR) field. We observe that the bulk RR field can give an effective 
potential which prevents the internal volume from collapsing. The flux gives a logarithmic 
bounce to the effective potential of the internal volume. 

In Sec. II, we consider the brane universe in the frame of Einstein gravity. We present 
a formalism for anisotropic expansion of brane universe and show that it is impossible to 
stabilize the internal volume with brane gases only. In Sec. Ill, we set up a supergravity 
formalism with brane gas and bulk RR field. For the anisotropic expansion, we consider the 
global rotational symmetry is broken down to SO(p) x SO(D — p) by (D — p)-brane gas. We 
decouple the equations of motion for each subvolumes of SO(p) and SO(D — p) and argue 
the possible stabilization. In Sec. IV, we conclude and discuss. 

II. ANISOTROPIC EXPANSION OF BRANE UNIVERSE 



In the previous work of us 11], we have shown that the extra dimensions should be 
wrapped with three- or higher- dimensional branes for the observed three large spatial di- 
mensions. However, it is not guaranteed that the compact dimensions can be stabilized with 
these branes. We start, after the thermal equilibrium of brane gas is broken, from the case 
when three dimensions are unwrapped and D — 3 dimensions are wrapped by gas of branes 
whose dimensions are less than or equal to D — 3. We assume that each type of brane gases 
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makes a comparable contribution to the energy momentum density. Then, we can study the 
evolution of a universe based on the following form of metric 

qmn = diag(-l, a 2 5ij, b 2 5 mn ), (1) 

where a 2 is the scale factor of the 3-dimensional space and b is the scale factor of the internal 
(D — 3)-dimensional subspace. 

The evolution of the universe with the above setting is equivalent to the asymmetric 
inflation in the brane world scenario [3] . We will consider a universe which can be described 
by the effective action 

S = J d ° +lx V9[^ - Vb+iGfo*)] , (2) 
where k 2 = l/Mf -1 , M* being the (D + l)-dimensional unification scale, and Vd+i(5'^) is 
the effective potential for the metric generated by the breakdown of the Poincare invariance. 
Considering the spatial section to be a D-dimensional torus T D , we write the metric as 

ds 2 = -n 2 dt 2 + J2 4(t)(dx k ) 2 . (3) 

k 

Technically, we retain the lapse n because we work on the action to reduce the action to a 
one-dimensional effective action. When we vary the reduced action with respect to the lapse, 
we obtain the Hamiltonian constraint. Later, we will set n = 1. Also we use the metric 
(J3]) instead of metric (pQ) to produce the correct combinational factors when the partial 
derivatives are taken over different scale factors. After all the derivatives are evaluated, we 
will take the identities a\ — ■ ■ ■ — a% — a and a^ = • •• = ajj = b. 

The effective potential Vd+i arising from the brane tensions can be written as 

D D 

y/gV D +i(9ixu) = n |a d JJ a k + T {p} JJ a k + JJ a k U (a k ) , (4) 

k=l {p} {p} k=l 

where is the tension of a D-dimensional space-filling brane which can be considered as 
a bulk cosmological constant, T{ p } is the tension of p-brane averaged over space, and U(a k ) 
parameterizes all contributions to the potential other than branes. For example, Fluxes, 
curvatures, and brane-brane interactions can give non-zero terms to U(a k ). We will show 
that U(a k ) plays an important role for the internal dimensions from collapsing to zero. In 
the summation over the brane tension, the notation {p} refers to the fact that there may be 
branes with the same p oriented along different directions. For example, 2-brane can wrap 
(12) cycle, (23) cycle et cetera. 



Now, we evaluate the Ricci scalar in Eq. ([3]) and reduce the problem to the motion of a 
particle in D dimensions. The one-dimensional effective action for the scale factors a k is 

/ ±J - 12k z I ' naj naj ) 

J k=l j=l 3 j=l ■> 

+ Ad + —^ X) T wII afe + C/ M' ( 5 ) 

lln=i a n {p} {p} 

where Vd = j d D x is the coordinate spatial volume, which contains the branes, and is 
dimensionless. In these units, the scale factors a k are dimension-full quantities, with dimen- 
sion M~ . If we vary this action with respect to n and a k , we obtain the following set of 
equations of motion after some straightforward algebra and setting n = 1 at the end, 

(E - E(^) 2 = 2^ 2 {Ad + U(a k )} + ET {P} l[a k , (6) 

j=l a 3 3=1 ^ i\.k=i a k {p} {p} 

!r + E^ = ^ ! T {A. + ^K)} + m E r «II°» 

a n ^a n a t D — 1 (D - 1) ]J k=1 a k ^ | p j 

K 2 d g 

+ m „ XIE^^-II * 

^ ^ -r-r 2 5f/ k 2 A df/ 

-pro }^ T {p} a n-^— \\ a k ~ k a n — + — — -> a k — . (7) 

We emphasize here that the partial derivatives are meant to be taken over different scale 



factors, even if some of the scale factors may coincide on the background. Note that the 

behaves like a counting operator acting on Yl{ P } a k- If 

' dai 



derivative operator cli-^ behaves like a counting operator acting on Yl{ P } a k- If a k 
depends on aj, the operator produces 1; otherwise, it gives 0. Hence, the operator Y2i 
simply counts the dimension of the brane. 

Now, we consider the special case when the global rotational symmetry is broken down 
to 5*0(3) x SO(D — 3) by effective (D — 3)-branes whose tension we denote Td_3. We ignore 
the stabilizing potential U(a k ) for simplicity. In this case, the equations of motion in Eqs. 
and (j7j) are simplified to, taking units in which 2/t 2 = 1, 

6( ^ + (B _ 3 )( B -4)(^ + 6p-3)^A D + ^i, (8, 
<i + 2(V + (Z)-3)^ = 7 ^ T + ^i T %i, (9) 



a V v 'ah D-l 2{D - 1) a 



- + (D — A)( h -Y + 3^ = ^ L_Zk±. (io) 

b { > K b> ab D-l 2(D-1) a 3 1 J 

The largest difference between the sizes a and b can be induced if a grows faster than b, 
as dictated by Eqs. and ffTUl) . The key parameters for controlling the relative rates of 
growth of a and b are their accelerations, not their velocities. To ensure that a exceeds b 
by many orders of magnitude, the sources in Eqs. (Q and (fTOl) must produce the slow-roll 
conditions for b, making its acceleration small or negative while keeping the acceleration of 
a positive. 

If the volumes of wrapped and unwrapped subspaces are defined as 

C = a 3 , ( = b D ~\ (11) 
the second derivative terms in Eqs. ([9]) and (TTOl) can be written as 

C - + ^ - 3 |^- + D ~ 2 Td - 3 } (12) 

l4^ D -^-2vh^y M 

The above equations of motion reduce to a simple form by removing the coupled first- 
derivative terms: 

£ D — 1 D ( 2{D — 1) C ' 1 J 

Since the left-hand side of Eq. (1141) is a function of £ while the right-hand side is a function 

of (, each side should be the same constant. Taking the constant as the parameter E, we 

have the following set of decoupled equations for ( and £: 

C - E+ 2(D-i) C ' ( 5) 

l = E + ^». (10) 

Note that the constant E is not an independent parameter. For the given configuration of 
branes, it is related to brane tension Td_3 and cosmo logical constant Ad. The constant E 
can be found by putting the solutions of (jl5]) and (JIB"]) into the original equation ffT2|) and 

To make the unwrapped subvolume ( expand forever, we require the constant E to be a 
positive definite. Then, the solution for the subspace ( is 

CHCo + C^-d, (17) 



where Co is the volume of unwrapped subspace at t — 0, 7 = y/~E, and C\ = 2(d-i) ^1T~- 
For given E, three types of solution for £ are possible depending on the bulk cosmo logical 
constant A D . First, for A# > 0, we have 



t = fa' nt , 7i = ^ + §^A D . (18) 

Since 71 > 7, £ expands faster than (. This means that the wrapped dimensions expand 
faster than the unwrapped dimensions. To determine the constant E, we put the solutions 
<Hm and (TTHD into (H 



2/1 , Ci\ , . Cis f D — 2 T D _ 3 



^ 1 + f )+Vn(1 + f )= HDfl + 2CD3i)Ti- <19) 
This gives the following relation between the constant £7 and other parameters in the unit 
of 2/t 2 = 1, 

7 2 + 77i = 3-^-, (20) 
(7 2 + 77i)Ci = |^y^-3- (21) 

There is no positive value of E which can suffice both of the above equations. This means 
that the solutions of the types ifPTj) and (118j) are not possible. 

Second, for A/? = — |Ad| < and E — ^Ef |Ad| > 0, the possible solution for £ is of the 
form 



t = fa'* t , ^ = ^-£-5^1. (22) 

In this case, £ expands slower than (. Repeating the same procedure as in the previous case, 
to check whether this solution can be realized, we have 

7 2 + 772 = -3 7 ^, (23) 

(7 2 + 772)^1 = |§^)^-3- (24) 

Since E > is a strict condition for the expansion of the unwrapped dimensions, it is obvious 
that the first condition cannot be satisfied. 

Finally, for A D = — \A D \ < and E - j^\A D \ < 0, we have 



4 = CneosM). ^ = x f -E + ^—^\A D \. (25) 
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In this case, the wrapped internal subvolume £ oscillates while the unwrapped volume ( 
expands monotonically. The solution seems quite similar to the one found by Patil and 
Brandenberger [3| in string gas cosmology in five dimensional spacetime without dilaton. 
They found that an oscillating type of solution, where the radius of extra dimension oscillates 
at the self-dual radius, is one of the possible trajectories. To check whether there is any 
constant E which gives this possibility, we substitute (fT7|) and (|25|) into (TT2]) 

^ + f ) + ^(1 + f ) ^3{dM + ^y^}. (26) 

It obvious that there is no nontrivial solution. We conclude, from the results for the three 
cases, that it is not possible to stabilize the internal dimensions with only branes. 

The above analysis with effective (D— 3)-branes is a special case of the dynamics described 
by Eq. (J7J). Other distributions of brane sources can give different possibilities. However, to 
have three large observed dimensions, we cannot neglect the potential U(ak) which represents 
the potential other than branes. Considering the motion of two subvolumes ( and £ as 
motion of particles in one-dimension effective potentials, the crucial point is that there is no 
bouncing potential for small £ for the third case in which we are interested. So the size of the 
internal subvolume is not bounded from below. This means that the internal subvolume can 
shrink to zero and can have even negative value, violating the assumption of the late-time 
brane gas approximation. 

III. BRANE GAS COSMOLOGY WITH BULK RR FIELD 

Having established the formalism and shown that it is not possible to stabilize the inter- 
nal volume with only brane gas, let us consider a model where the potential U{a>k) comes 
from bulk flux. The model we will consider is type II string theory compactified on T 9 . Also 
we consider the late stage of the BV scenario where the radii and curvature scales are grown 
larger than the ten-dimensional Planck length. Since winding numbers are topologically 
conserved, only branes with opposite orientations can cause unwinding. This unwinding 
interactions are very short ranged since this happens when the world volumes of branes 
and antibranes physically intersect as in BV scenario When the radii of the unwound 
directions have grown enough, the brane and antibrane densities will be diluted so that they 
cannot find each other to intersect. Then we can neglect the brane-antibrane annihilation 
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since this term will be turned off as the transverse dimensions expand 0, Q, Super- 
gravity is a good approximation with the growing radii and falling temperature. In brane 
gas cosmology it is assumed that brane mode fluctuations are small and only winding modes 
dominate the cosmological evolution. So we consider the wrapped D-branes, which do not 
interact with each other, for the matter part of the Einstein equation. This corresponds 
to a classical approximation of the brane dynamics. Assuming that the brane gases are 
homogeneous, we can take the average of the contributions of all kinds of branes to the 
energy-momentum tensor. 

In the point of string theory, the gravitational interaction is described by the coupled 
system of the metric and dilaton. Dilaton plays an important role in the large-small sym- 
metry of string theory called T-duality. The initial big-bang singularity can be resolved in 



terms of T-duality [l, 2]. We start from the following bulk effective action of type II string 
theory 

= / d D+1 x^-g[e-^{R + 4(V0) 2 } - ±Ff - V(<f>)] , (27) 

where is the dilaton field, Fi is an RR (Z)-form field strength and V(<p) is the dilaton 
potential. The dilaton potential is a product of nonperturbative effects. For simplicity, we 
consider only flat potential which equivalently is a cosmological constant. 

The matter contribution of a single brane to the action is represented by the Dirac-Born- 
Infeld (DBI) action of p-brane 

S P = -T p J dF+^e-tyJ -det(^ + B, u + MF~), (28) 

where is the induced metric to the brane 

dx M dx N 

9nu = 9mn q ^ ■ (29) 

Here M, N are the indices of bulk spacetime and fi, v are those of brane. is the induced 
antisymmetric tensor field and F^ u is the field strength tensor of gauge fields living on 
the brane. The fluctuations of the brane coordinates and other fields within the brane are 
negligible when the temperatures are low enough. So we neglect B^ v and F^ v terms below. 
The total action can be written as 

S to t = S h + J2S{ p} - (30) 

M 
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Since we are interested in the late time cosmology, we have to work in the Einstein frame 
defined by gf^ N = G D - 1 ^'gf tN . The formalism and physical picture for the investigation 
of metric perturbations have been developed in the Einstein frame where one can use the 
standard simple Einstein equations. Since the weak coupling limit (i.e. g s = « 1) is one 
of the basic assumptions of string/brane gas cosmology, we ignore the running of dilaton 
assuming that dilaton can be stabilized. When the dilaton is not running, the two frames 
are equivalent since the contribution of constant dilaton can be absorbed in the redefinition 
of brane tension and RR field. Then, the action can be written as 



S% = ^2 / d D+1 x^\R - ±- } F? - A D ] , (31) 
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S P = -T P J d p+1 ^-detgz u . (32) 

We drop the superscript E from now on. 

In the point of bulk theory, the energy momentum tensor of a single Dp-brane has a delta 
function singularity at the position of the brane along the transverse directions 

^~ gT MN = ^ J d P+^^gg^ d ^ X M duX N 5[x _ ^ (33) 

In cosmological setting, it seems natural to take a gas of such branes in a continuum ap- 
jroximation and this smooths the singularity by integrating over the transverse dimensions 

The induced metric on the p-brane follows from its embedding in the bulk. We choose 
the static gauge £ M = x M (/i = 0, 1, . . . ,p) and suppose the simple embedding x % = x l {t) for 
i = p + 1, . . . , D. The induced metric on the brane is, 

= diag{-(l - v 2 ), a 2 , a%,..., aj}, (34) 

where v 2 = J2^ =p+l x % ii > 0. Note that the brane time, defined by dr = dt\/l — v 2 , is not 
the same as the bulk time if the brane is moving. For the RR field we take the ansatz 

Fn u -,Ni = nVle Nl> ... tNl _ 1 'VoA{iy(t), (35) 

where N±,--- , iV;_i ^ 0. Then the Bianchi identity 

V[nF Ni ,.. m = 0, (36) 
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is automatically satisfied since Ai(t) is a function of t only. 

If we consider the static brane, i.e., consider the brane does not move in the transverse 
direction (v = 0), the brane time is the same as the bulk time. Then the one-dimensional 
effective action for the scale factors can be written as 

D D D 



• ; k=l j=l J j=l 3 

~ Ad + fr°~ — E T {piH Gfc + i E 

ll„=i a n f v x M m 1 1.7=1 a j 



(37) 



to to W 1^=1^ 

We take variation of this action with respect to n, aj~, and A^\. After some straightforward 
algebra setting n — 1 and 2k 2 = 1 at the end, we obtain the following set of equations of 
motion: 

D . D . , ^ A 2 



v ' llfe=l « m=l {p} {p} m=l {£} 1 lj=l J 



"^{fP^T E T w a ™^- 11^ + E fpr^}> ( 39 ) 

Z L llfc=l«re {p} n {p} ^ {1} llj=i a 3 



d Aj 



dt 



{n-n^} 

fe=l 1 lj=l "j 



0. (40) 



Note that, comparing Eqs. ( 138]) and ( l39l) with Eqs. ([6]) and (171), the potential £7 comes from 
the Ay\ term originated from bulk RR field. 
The solution of (|40p is given by 



J- 1 „2 



^{i} rij=i 



CI': 



Am = ^ 3=1 \ (4i) 

llk=l fl re 

with an integration constant Q{i}. Substituting (14T1) . the equations of motion ff38]) - (|39]) can 
be written as 

i=i a J ,-=i a ^ llfc=i°fc {p } { p } z w llfc=i a K 



n 



^ u y > K ih=i a k {p} {p} {l} iU=iV 

-5feg r » a *^B a * + ?-^5 0!o ^ } ' (43) 

For the anisotropic expansion of dimensions, we consider the global rotational symmetry 
is broken down to SO(p) x SO(D — p) by (D — p)-brane gas. What we are interested is the 
case with p = 3. But we consider the general case to see the dependence of dimensionality 
on branes and fluxes. For simplicity, we introduce two kinds of bulk RR potential and brane 
which are characterized by Q p , Qd- p , T p and Td_ p . Taking the scale factor for p subspace 
as a and D — p subspace as b, the equations of motion are given by 

p(p - 1)(-) 2 + (D — p)(D — p — 1) A 2 + 2p(D -p)- 
a b ao 

= Ad + + lEzIL + 1 Q p + I^z£ (44) 
a , . / n ^ Ad D-p-2 T p D - p + 1 T D _ P 

- + (p-i)(-) +^-^^ = ^1-^317^ + ^^17^ 
»(D- P + i) + ig^_ p (p-p) 2 -i q 2 

2(D-1) a 2 P 2(D-1) 6 2 ( d -p)' 1 j 

iU(p-p-i)A 2 + P ^ = -^- + p+1 p ~ 2 Td ~ p 

b { P A b> P ab D - 1 2(D - 1) b D -P 2(D-1) a? 

(D-p)(p-i) + i Ql _ p 2 -i Ql- P ( v 

2(D-1) 6 2 ( d -p) 2(D-1) a 2 f ' 1 J 

If we define the two volumes of SO(p) and SO(D — p) as 

C = a?, £ = b D ~ p , (47) 

the second derivative terms in Eqs. (T4"5!) and (jlBjl can be written as 
C a f A D D- P +1Td_ p D — p — 2T p 



c a ^d-i 2(d-i) c 2(d-i) e 

P (D-p-i) + ig^_ p (d-p) 2 -i q 2 

2(£>-l) C 2 2(D-1) e 1 
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£ , tt_ {D p) f Ad , P+l T p p-2 T D _ P 



e a ' "l£>-l 2(£>-l)£ 2(D-1) C 

(D-p)(p-l) + i g; _ qi_ a 
2{D-l) i 2 2(D-1) C 2 J' 
Decoupling the variables by introducing a constant as in the previous section, we have 
C = E p(D-p+l) + (D-p)(p-2)T D _ p 

C 2(D-1) C 

p 2 (jD _ 1} + + (jD _ p)(p2 _ 1} 

+ 20D^1) C 5 "' (50) 

I _ „ , £~2p , (£>-p)(p+l)+p(D-p-2) T p 

e" E + ^r AD + 2p^i) j 

(g - p)2( p _ 1) + g _ g + g((g _ p) 2 1) Q2 

+ 2(D-1) e 2 ' 1 ' 

As a consistency check, one can recover Eqs (fT5T) and (TT6l) by taking p = 3 and T p = Qd- p = 

Q p = 0. Note that the constant E is not an independent parameter. For given configuration 

of branes and RR fields, it is related to A D , T p 's and Q p 's by putting the solution obtained 

from the decoupled equations (l50i) and ( 15T1) into the original ones (jlSjl and ( 1491 ). 

Since Eqs ( 1501) and (15TT) are non-linear because of the Q 2 terms, finding a closed form 

of analytic solution seems very difficult. So we search the possibilities of the solutions. 

We analyze the behavior of the two subvolumes by considering the effective potential £ = 

— dVeB ^ £ = — dV e«{£) _ xhe effective potentials are obtained as 



where 



Kff(C) = -|C 2 -CiC-c 2 lnC, (52) 
V eS (0 = -^e-di£-d 2 ln£, (53) 



c = E, 

= P( D- P+1) + ( D -P)(P~2) 

2(D - 1) p y ' 

p\D-p-l)+p+ (D-p)(p 2 - 1) ^ 2 

D ~2p k 

do = E+ — _ ^ A g , 

2(£) - 1) p V 7 

(g - p) 2 (p - 1) + g - p + p((£) - p) 2 - 1) ^ 2 
d2 - 2{D~l) Qp ' 
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FIG. 1: Typical shape of effective potential V e g(() for unwrapped subvolume for E > 0. The plot 
is for cq = 1, ci = 1, and c 2 = 0. 



Now let us examine the possible combinations of brane gas and background RR flux. To 
see the contribution of RR potential to the evolution of brane universe, consider the simple 
case when Td_ v and Q p are nonzero with p = 3. This corresponds to the case that we 
considered in Sec. II with the inclusion of 4-form RR flux in the unwrapped dimensions. 
For effective potential of the observed subvolume V e s((), we have 

c = E, 

AD -9 

ci = ap^if"-*' (56) 
c 2 = 0. 

The shape of the effective potential for E > is given in Fig. 1. In this case the equation of 
motion for ( (Eq. (150]) ) is the same as Eq. ffT5l) . So ( expands forever as far as E is positive. 
For effective potential of the internal subvolume V^g (£), we have 

r, D - Q ^ 
d = E + jj—^^d, 

d l = 0, (57) 
_ 5(£>-3) 2 + £>-6 ^ 2 

We need the condition = — |Ad| < and E — |Az?| < as in Sec. II to have the 
confining behavior for large £. To have a bouncing potential for small £, we need c?2 > 
and this condition is automatically is satisfied. The shape of the effective potential for this 
case is given in Fig. 2. The wrapped subvolume £ will oscillate around the minimum of the 



effective potential £ min = ^/—d 2 /d = ^/d 2 /\d \. 
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FIG. 2: Typical shape of effective potential V e s(£) for wrapped subvolume for Ad = — \Ad\ < 
and E - %^\A D \ < 0. The plot is for d = -1, d 1 = 0, and d 2 = 1. 



\ 




FIG. 3: Typical shape of effective potential V c fr(C) f° r unwrapped subvolume when there exists 
(D — 3 + l)-form RR field in the internal dimensions. The plot is for Co = 1, c\ = 1, and C2 = 1 

The existence of (3+l)-form RR field induces a logarithmic bounce on the effective po- 
tential of the subvolume perpendicular to it. This bounce prevents the internal subvolume 
from collapsing to zero size. Similarly, the existence of (D — 3 + l)-form RR field in the 
internal dimensions can induce nonzero logarithmic bounce for (. However, as far as the 
condition E > is satisfied, the overall behavior of ( is essentially the same (see Fig. 3). 



IV. CONCLUSION 

In the framework of brane gas cosmology, we considered the asymmetric evolution of the 
universe. We showed that, provided the dilaton is stabilized by some mechanism, brane 
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gas and RR potential can stabilize the internal dimensions. For a combination of (D — 3)- 
brane gas wrapping the extra dimension and 4-form RR flux in the unwrapped dimensions, 
it is possible that the wrapped internal subvolume has an oscillating solution around the 
minimum of the effective potential while the unwrapped subvolume expands monotonically. 

One key point of this paper is bulk RR field. The role of this field is to prevent the 
internal subvolume from collapsing to zero. It gives a logarithmic bounce to the effective 
potential for small values of internal subvolume. Another is the decoupling of the observed 
and internal dimensions by introducing an integration parameter. With this decoupling, we 
could treat the problem analytically in terms of two subvolumes rather than scale factors. 

In our analysis we considered only brane gases for the matter part of the Einstein equation. 
However, there can be supergravity particles corresponding to bosonic and fermionic degrees 
of freedom. One can ignore the massive modes since these will decay quickly. For the 
massless supergravity particles, if we assume they are homogeneous and isotropic, we can 
take a perfect fluid form of energy momentum tensor 

T*f = diag(-ps,ps,-.. ,Ps). (58) 

It is well known that this term always tends to drive uniform expansion. So the main result 
of our analysis will not be changed even if we include the massless supergravity particles. 



171 ]. where the extra dimensions are 



To have more complete stabilization like the one in 
stabilized at the self-dual radius, one can include the effect of other sources. For example, 
including the effect of momentum modes or running dilaton may give damped oscillatory 
solution in which £ is driven to the minimum of the potential. 

We would like to point out that the configuration in this work is just one of the possibilities 
to achieve the anisotropic expansion. To make anisotropic expansion between the scale 
factors of the observed 50(3) and those of SO(D — 3), we considered the configuration 
where a gas of positive-tension (D — 3)-branes wrap the internal dimension. An alternative 
way to achieve the same anisotropic expansion is to use a gas of negative-tension 3-branes 
wrapping the observed dimensions. Other combinations of brane gas and RR flux may give 
more desirable result. 
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